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I.  Introduction 

This  project,  entitled  “Computations  of  Optimal  Controls  and  Designs  for  Distributed 
systems  in  Optics  and  Elasticity”,  was  funded  by  two  separate  grants  from  AFOSR  (Grant 
#88-0091)  and  NSF  (Grant  #DMS-871S510).  AFOSR  has  provided  additional  release 
time  for  Drs.  G.  Chen  and  J.  Zhou,  the  two  Pi’s  for  this  project,  while  NSF  has  provided 
fund  for  the  acquisition  of  a  SUN  3/60  Workstation. 

During  the  support  period  January  through  December  1989,  the  Pi’s  have  studied 
various  questions  of  controls,  designs  and  computations  using  the  boundary  element  method 
(BEM)  for  distributed  systems  governed  by  the  following  equations  in  optics  and  elasticity: 

1)  the  Helmholtz  equation  on  an  exterior  domain; 

2)  the  eigenvalue  problem  for  the  Laplacian; 

3)  the  biharmonic  elastostatic  Timoshenko  plate  equation;  and 

4)  three  dimensional  solid  mechanics  equations. 

Numerical  software  has  been  accumulated  and  computer  graphics  has  been  successfully 
developed.  So  far,  two  papers  have  been  published.  Three  others  are  being  prepared  for 
publication.  See  §11  and  §111  below. 

This  research  project  has  been  consolidated  with  our  AFOSR  Grant  87-0334  and  will 
be  carried  on  further  in  that  grant. 

§11.  List  of  Publications 

The  following  research  papers  have  been  published  under  the  auspices  of  this  grant: 

[1]  Computing  optimal  boundary  controls  of  a  plate  by  the  boundary  element 

method  (authors:  G.  Chen  and  J.  Zhou) 

Status:  presented  at  the  26th  IEEE  Control  and  Decision  Conference,  Los  Angeles, 

December  1987,  and  appeared  in  the  Proceedings,  pp.  992-996. 

To  our  knowledge,  this  is  the  first  successful  computation  of  the  biharmonic  Timo¬ 
shenko  plate  equation  by  the  boundary  element  method.  Optimal  boundary  controls  in 
terms  of  displacement  and  angle  are  also  computed. 
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[2]  Minimizing  the  reflection  of  waves  by  surface  impedance  using  boundary 
elements  and  global  optimization  (authors:  G.  Chen,  T.J.  Bridges  and  J.  Zhou) 
Status:  published  in  Wave  Motion  10  (1988),  pp.  239-255.  The  material  was  also 
presented  at  an  unscheduled  talk  in  the  Session  on  Control  of  Distributed  Parameter 
Systems  in  the  26th  IEEE-CDC  in  Los  Angeles  in  December  1987. 

In  this  paper,  we  used  the  boundary  element  method  and  global  optimization  tech¬ 
niques  to  study  the  “stealth  problem",  i.e.,  the  minimization  of  reflected  radar  or  sonar 
waves.  All  of  the  computer  codes,  including  the  Helmholtz  equation  solver  on  an  exterior 
domain  and  the  global  optimization  algorithms,  were  developed  by  us.  Numerical  results 
have  been  benchmarked  and  have  conformed  to  the  physical  circumstances. 

§111.  Papers  in  Preparation 

The  following  papers  are  presently  being  written  for  submission  to  scientific  journals. 
Preprints  will  be  sent  to  the  Program  Manager  at  AFOSR  once  they  become  available. 

[1]  Boundary  integral  equations  and  boundary  element  computations  of  the 
Timoshenko  plate  equation  with  applications  to  boundary  shape  control 
This  is  a  long  technical  paper  where  we  have  formulated  the  boundary  integral  equa¬ 
tions  for  the  elastostatic  biharmonic  plate  equation  subject  to  four  types  of  combinations 
of  displacement,  angle,  bending  moment  and  shear  boundary  conditions.  Boundary  ele¬ 
ment  computations  have  been  successfully  carried  out  and  tested.  The  methods  are  then 
applied  to  boundary  shape  controls  and  numerical  modelling  of  the  control  system. 

We  expect  this  work  to  be  complete  by  early  May  1989. 

[2]  The  boundary  element  method  for  shape  control  of  distributed  parameter 
systems 

Drs.  G.  Chen  and  J.  Zhou  have  been  invited  by  Professor  John  L.  Junkins  to  write  a 
survey  chapter  for  a  new  volume  entitled  “ Mechanics  and  Control  of  Space  Structure"  in 
the  AIAA  “Progress  in  Astronautics  and  Aeronautics  Series",  for  which  Professor  Junkins 
is  the  editor.  This  volume  will  summarize  the  state  of  arts  of  active  control  of  space 
structure. 
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The  due  date  of  the  first  draft  of  the  manuscript  is  June  1,  19S9. 

[3]  The  eigenvalue  problem  and  eigenshapes  for  the  Laplacian 

Traditionally,  the  computation  of  eigenvalues  and  eigenfunctions  of  the  Laplacian 
on  nonrect angular  domains  is  done  by  the  finite  element  method,  which  involves  a  large 
amount  of  coding  work,  numerical  quadrature  and  computer  storage.  We  have  developed 
boundary  element  methods  and  algorithms  for  studying  such  problems,  which  are  accurate 
and  fast.  Computer  graphics  for  eigenshapes  has  also  been  successfully  produced,  which 
displays  remarkable  accuracy  in  comparison  with  the  benchmarks. 

This  work  is  expected  to  be  complete  by  October  1989. 

[4]  The  boundary  element  method  and  shape  control  for  3-D  solid  mechanics 
The  SUN  3/60  Workstation  acquired  with  NSF  Grant  DMS  87-18510  has  been  in¬ 
stalled  in  our  office  since  May  1988.  We  are  now  able  to  run  the  machine  with  UNIX  and 
have  become  familiar  with  SUN’s  3-D  color  graphics  package  for  our  visualization  purpose. 
For  3-D  solid  mechanics,  the  governing  equations  are 

3  o 

T  o — 0ij(x)  =  /«(*)»  x  6  fi  C  R3,;  =  1,2,3 

UdZi 

where 


<7,j(x)  =  the  linear  stress  tensor, 

=  A[V  •  u(x)]6ij  +  2 pejj(z), 

£ij(x)  =  the  linear  strain  tensor, 
u(x)  =  the  displacement  vector, 

A ,  fi  are  Lame’s  constants, 

/,(.r),  i  =  1,2,3,  are  components  of  body  force. 

A  boundary  element  package  has  been  developed.  To  illustrate  our  progress  and  achieve¬ 
ment,  we  include  the  graphics  from  a  simple  numerical  example,  using 
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/*(/l,U)3  0 


n  =  the  unit  ball 
r  =  the  boundary  traction  force 
=  (ti,t2}  r3) 

rj  =  r2  =  0  on  dCl 

{—(1  —  sinci),  0  <  <t>  <  rr/2 

,  0  is  the  latitudinal  angle. 

1  —  sin  0,  j  <  o  <  t 

Therefore  the  ball  is  subject  to  a  vertical  compression  force.  In  Fig.  1  (a)  and  (b),  we  see 
the  unit  ball,  whose  surface  is  discretized  into  30  x  40  =  1, 200  piecewise  constant  boundary 
element.  In  Fig.  2  (a)  and  (b),  we  see  that  the  bail  has  been  compressed  into  an  ellipsoidal 
shape.  In  Fig.  3  (a)  and  (b),  we  see  that  the  ball  has  been  further  compressed  into  an 
“English  muffin”  if  the  vertical  compression  force  (*)  is  tripled. 

This  research  of  boundary  element  computations  and  applications  to  solid  mechanics 
control  was  presented  by  G.  Chen,  L.  Ji  and  J.  Zhou  at  am  Invited  Session  on  Inverse 
Problems  at  the  27th  IEEE-CDC,  Dec.  7-9,  1988,  in  Austin,  Texas. 

A  research  paper  is  being  prepared  for  publication. 
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(b) 


Fig.  I 


In  (a)  and  (b),  the  surface 
30  x  40  piecewise  constanc 


of  the  unit  sphere 
boundary  elements 


is  discretized 


into 


subject  to  a  compressive  traction  force 


wa  ->1  3  !  "(l~sin  0<  b<TT/ 2  ,  ,  ,  , 

m*'0”-'  •  ,  0  =  longitutinaL  angle,  4>=latitudinal  angli 

l-sin  ;  ,  t/2<0<^ 


3  The  deformed  unit  bail  takes  the  shape  of  a  "muffin",  where  the  surface 
is  subject  to  a  compressive  traction  force  three  times  of  that  in  Fig.  2. 
Here  the  extent  of  deformation  actually  has  exceeded  the  linear  elasticity 
limit.  The  purpose  of  this  figure  is  to  illustrate  the  dramatic  effects. 


§IV.  Personnel 

Dr.  G.  Chen,  the  PI,  supervises  the  overall  conduct  of  the  grant.  Dr.  J.  Zhou  is  the 
co-PI.  The  following  people  are  collaborators  associated  with  the  research  effort: 

T.J.  Bridges  -  Assistant  Professor.  Department  of  Mathematics,  Worcester  Poly¬ 
technic  Institute,  Worcester,  MA  01609 

L.  Ji  -  Ph.D.  student,  Department  of  Mathematics,  Texas  A&M  University,  Col¬ 
lege  Station,  TX  77843 

§V.  Activities 

Drs.  Chen  and  Zhou  have  given  the  following  talks  on  their  collaborated  work: 

[1]  The  26th  IEEE-CDC,  Los  Angeles,  CA,  Dec.  1987. 

Subject  matter:  Computing  optimal  boundary  controls  of  a  plate  by  the  boundary 
element  method. 

Individual  involved:  Prof.  J.E.  Lagnese,  Georgetown  University  and  Prof.  M.C. 
Delfour,  Univeristy  of  Montreal,  Canada. 

[2]  Workshop  on  Numerical  and  Experimental  Aspects  of  Distributed  Parameter  Control, 
Center  for  Mathematical  Sciences,  University  of  Wisconsin,  Madison,  WI,  May  1988. 
Subject  matter:  The  boundary  element  method  for  computing  boundary  controls  in 
optics  and  elasticity. 

G.  Chen  is  a  major  invited  speaker  in  the  Workshop,  wherein  J.  Zhou  gave  a  demon¬ 
stration  of  coirlputer  graphics. 

Individual  involved:  Prof.  D.L.  Russell,  University  of  Wisconsin. 

[3]  Conference  in  honor  of  the  65th  birthday  of  Regent  Professor  Lawrence  Markus,  Dept, 
of  Math.,  University  of  Minnesota,  Minneapolis,  April  1988. 

Subject  Matter:  Boundary  control  of  a  thin  elastostatic  plate  by  the  boundary  element 
method. 

G.  Chen  was  an  invited  keynote  speaker  for  the  conference. 

Individual  invoked:  Prof.  E.B.  Lee,  Electrical  Engineering  Dept.,  University  of 
Minnesota. 

[4]  The  27th  IEEE-CDC.  Austin.  Texas.  Dec.  19SS. 
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Subject  matter:  The  boundary  element  method  for  shape  control  in  3-D  solid 
mechanca. 

Individual  involved:  Prof.  G.  Crosta,  University  of  Milan,  Italy. 


